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ABSTRACT 



We present a self-consistent three-dimensional Monte-Carlo radiative transfer model of the stellar and dust emission in 
the Milky- Way, and have computed synthetic observations of the 3.6 to 100 (im emission in the Galactic mid-plane. In 
order to compare the model to observations, we use the GLIMPSE, MIPSGAL, and IRAS surveys to construct total 
emission spectra, as well as longitude and latitude profiles for the emission. The distribution of stars and dust is taken 
from the SKY model, and the dust emissivities includes an approximation of the emission from polycyclic aromatic 
hydrocarbons in addition to thermal emission. The model emission is in broad agreement with the observations, but a 
few modifications are needed to obtain a good fit. Firstly, by adjusting the model to include two major and two minor 
spiral arms rather than four equal spiral arms, the fit to the longitude profiles for \£\ > 30° can be improved. Secondly, 
introducing a deficit in the dust distribution in the inner Galaxy results in a better fit to the shape of the IRAS longitude 
profiles at 60 and 100 (jm. With these modifications, the model fits the observed profiles well, although it systematically 
under-estimates the 5.8 and 8.0 /im fluxes. One way to resolve this discrepancy is to increase the abundance of PAH 
molecules by 50% compared to the original model, although we note that changes to the dust distribution or radiation 
field may provide alternative solutions. Finally, we use the model to quantify which stellar populations contribute the 
most to the heating of different dust types, and which stellar populations and dust types contribute the most to the 
emission at different wavelengths. 

Key words. Radiative transfer - Galaxy: structure — Infrared: ISM 



1. Introduction 

Although it has been common to characterize the large- 
scale Galactic diffuse emission by treating each component 
of the Galaxy individually - the neutral gas distribution, 
HII region maps, the magnetic field structure - this emis- 
sion is usually caused by a convolution of different combina- 
tions of these components. Far-infrared dust emission, for 
example, traces the dust distribution heated by stars; radio 
free-free emission depends specifically on the local density 
of main-sequence massive stars; radio synchrotron emission 
traces the local cosmic ray electron density (and energy 
spectrum) convolved with the magnetic field structure; and 
so on. Since our knowledge of the distribution of any one 
of these components is still rather incomplete, modeling 
emission that arises from the coupling of two or more of 
them is a daunting task, requiring numerous assumptions 
to develop a complete model that can be compared with 
observations. 

The value of such an exercise lies in identifying the fac- 
tors that have the largest influence on determining the 
properties of the observed emission, with the hope that 
eventually all of the different diffuse emission measurements 
can be combined to yield a more complete picture of the 



Galaxy. In this paper, we present models of the stellar and 
diffuse infrared emission in the Galactic plane. Although in- 
vestigations of the diffuse infrared emission have been car- 
ried out in the past, to our knowledge this is the first work 
to focus on the mid-infrared (3.6-8.0 /im) diffuse emission 
from polycyclic aromatic hydrocarbons (PAHs), which is 
much more sensitive to the distribution of the young stars 
than the longer wavelength dust emission probed by IRAS 
and COBE/DIRBE. 



The large-scale distribution of stars and diffuse dust 
emission in the Milky- Way has been the subject of numer- 
ous studies. On one hand, powerful three-dimensional stel- 
lar distribution models have been developed to simulate the 
inventory of stellar populations as a function of location in 
the Galaxy and position on the sky. The Besanqon model 
(iRobin fc Crezd (l986t iBienavme et~all 119871: IRobin et al l 
| 1996l I2003l)7~the SKY model (|Wainscoat et all 119921: 



Cohed Il993l jilli [111), and the TRILEGAL model 



(jOirardi et al.1 2005) are three of the most widely known 
such models aimed at reproducing the stellar populations 
of the Galaxy. Such models typically include various compo- 
nents to describe the Galaxy, such as a bulge, halo, disk(s), 
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spiral arms, and so on, and include stars at various masses, 
metallicities, and stages of evolution. 

On the other hand, analytic models of the diffuse dust 
emission in the Galaxy have also been developed to re- 
produce t he large-scale mid- an d far-infrared emission. For 
example, ISodroski et all (|1997t ) used COBE/DIRBE data 
from 12— 240 /mi to determine the best-fitting abundances 
and temperature of large dust grains, the abundance of 
very small transiently heated dust grains and PAHs, and 
the energy density of the interstellar radiati on field as a 
functi on of Galactocentric radius. Similarly, iDavies et al.l 
(|1997t) fit a model for the emission from cool (18-22K) 
dus t to the COBE/DIRBE far-in frared data. More recently, 
the iPlanck Collaboration! (|201lf) fit the observations of the 
Galaxy by t he Planck satel lite, using a model similar to 
that used by ISodroski et al.l 

IDrimmel (|2000l ) carried out a qualitative comparison 
of COBE/DIRBE K-band and 240 /jm observations, and 
found that while the near-infrared shows evidence only for 
two stellar arms, the far-infrared observations are consis- 
tent with four arms (in agreement with the influential HII 
region study of lGeorgelin fc Georgelin|[l970D , leading to the 
conclusion that the main Galactic potential is two-armed, 
but that the structure of the dust and gas responding to the 
potential is more com plex, and can be ad e quate ly charac- 
terized as four-armed. IDrimmel fc Spergel (|200lD extended 
this by developing a quantitative model to fit the above 
data, which consisted of a parametric model of the dust 
distribution and associated fa r-infrared emiss i on. T his dust 
distribution was then used bv IDrimmel et al.l (|2003l ) to pre- 
dict the three-dimensional K-band extinction. 

What has been lacking in previous models is a self- 
consistent treatment of the source of the dust heating - 
the photons produced by stars of various temperatures and 
luminosities - and the target of the heating - grains of dif- 
ferent sizes and PAHs - distributed throughout the disk of 
the Galaxy. In this paper we present a radiative transfer 
model of the Galaxy at infrared wavelengths, considering 
both stellar and dust emission. The model is self-consistent 
in the sense that the heating of the dust does not follow 
an analytical prescription, but is instead directly computed 
from the stellar populations. The main aims of this paper 
with this model are the following: 



— To determine whether an existing Galactic model for 
stellar populations and dust - here the SKY model - 
us ed in conjunc t ion w ith the dust properties determined 
bv lDraine fc Lil (|2007f ). which include transiently heated 
very small grains and PAH molecules, adequately repro- 
duces the stellar and diffuse emission seen at mid- and 
far- infrared wavelengths. 

— To determine whether any modifications are needed to 
this model to reproduce the observations, and whether 
any new insights on Galactic structure can be gained. 
However, we stress that the aim is not to provide an 
exhaustive parameter space study, as given the number 
of stellar populations and components, the number of 
parameters describing the model is of the order of sev- 
eral hundred. Rather, the aim is to determine possible 
and realistic modifications that improve the model. 

— To determine the relative contributions of resolved and 
unresolved stellar flux, emission from PAH molecules 
and dust grains, and scattering at infrared wavelengths. 



— To determine the relative importance of various stellar 
populations in heating dust grains and exciting PAH 
molecules. 

In Section we describe the observational datasets 
used in this paper. In we describe the initial model used, 
the dust properties, and the radiative transfer code. In <21 
we show results for both the initial model f ^4.1[) as well as 
an improved model ( tj4.2p . In |J5]we examine the main con- 
tributors to the observed flux and to the heating of the dust 
(SSHJ), we study how much IRAC stellar flux is likely to be 
unresolved f H5.2[) . and we generate images of the Galaxy 
from an external viewpoint ( tj5.3[) . Finally, in fJ5] we sum- 
marize our findings. 

2. Observations 

In this paper w e mak e use of the Spitz er GLIMPSE 



(jBeniamin et al 



_ 12003 IChurchwell et aD l2009h and 

MIPS GAL (jCarev et al.l l2009D surveys of the Galactic mid- 
plane, specifically GLIMPSE/MIPSGAL I, which cover 
10° < \£\ < 65° and |6| < 1°, and GLIMPSE/MIPSGAL II, 
which fill in the region for \£\ < 10°, with |6| < 1° for 
\£\ > 5°, |6| < 1.5° for 2° < \£\ < 5°, and |6| < 2° 
for \£\ < 2°. The total area covered by these surveys is 
274 deg 2 . 

The GLIMPSE surveys wer e carried out us ing the 
InfraRed Array Camera (IRAC; iFazio et al.ll2004l) at 3.6, 
4.5, 5.8, and 8.0 fim. The MIPSGAL surveys were carried 
out usin g the Multiband Imaging Photometer for Spitzer 
(MIPS; iRieke et al.l I2004D at 24 and 70 /jm, though only 
the processed 24 /im data are available and are included in 
this study. 

In addition to the Spitzer data, we include data from 
the Improved Reprocessing of th e IRAS Survey (IRIS; 
iMiville-Deschenes fc Lagachd 120051) at 60 and 100 /im to 
provide further constraints on the dust emission. 

At the different wavelengths considered, the total emis- 
sion in a given direction consists of stellar emission re- 
solved into point sources, unresolved stellar emission, and 
diffuse interstellar emission. The relative contribution of 
these three types of emission change both as a function 
of wavelength (diffuse emission increases with wavelength), 
and position in the Galaxy (unresolved stellar emission de- 
creases with increasing angle from Galactic center). Here, 
longitude and latitude profiles were constructed for the four 
IRAC bands, the MIPS 24 /im data, and the IRIS data, by 
resampling the full data into 3' by 3' bins, including the re- 
solved stellar flux. The re-binned data were then collapsed 
into selected longitude and latitude profiles to compare the 
observations to the radiative transfer model presented in 
21 In this paper, we only use the rectangular region for 
\£\ < 65° and |b| < 1° to facilitate computing longitude and 
latitude profiles. 



3. Model 

Our self-consistent radiative transfer model of the stellar 
and dust emission at infrared wavelengths consists of two 
key ingredients: a model for the distribution of the stel- 
lar populations, and a model for the characteristics and 
distribution of the dust. The initial distributions used are 
described in fcl3.ll and in fc|3.2l but these are then modified 
in i j4.2l to provide a better fit to the observations. The dust 
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properties assumed are described in 
diative transfer model is described in 



Finally, the ra- Table 1. Spiral Arm Parameters 



3.1. Stellar populations 

To model the distribution of variou s stellar populations, we 
implemented the S KY model from Wainscoat et al.l (|1992t 
hereafter W92) and ICohenl I1994I . Il995h . 7 rhe choice 

of the SKY model over more recent models such as the 
Besancon or TRILEGAL models is that the former pro- 
vides a prescription for spiral arms, which are essential for 
the present work, since we expect the diffuse PAH emis- 
sion at IRAC wavelengths to depend on the distribution of 
massive stars, which are concentrated in the spiral arms. In 
addition, the three-dimensional distribution of the stellar 
populations in the SKY model is separated into spectral 
classes, making it straightforward to represent in a radia- 
tive transfer model. 

The SKY model includes five components to describe 
the Galaxy: an exponential disk, a bulge, a halo, spiral 
arms, and a ring. We do not include the halo: as it is defined 
in W92 by a projected rather than a spatial distribution, it 
is more difficult to include in the radiative transfer model. 
However, it has a very low n ormalization factor relative to 
the disk (l:1250: [CoheJl995l ) and would have no noticeable 
impact on the Spitzer or IRAS observations in the Galactic 
mid-plane. The remaining four main components are briefly 
described in the following sections. 

The model includes 87 different types of Galactic 
sources, including pre-main-sequence stars, main-sequence 
stars, giants, asymptotic giant branch (AGB) stars, super- 
giants, and planetary nebulae. Each component includes a 
normalization factor (pd, Pb, Pa, Pr), and in some cases 
other parameters that depend on the source type, denoted 
S. For example, O, B, and A stars are absent from the bulge 
(typically composed of more evolved stars), while lower- 
mass main sequence stars are absent from the spiral arms 
(typically composed of relatively young stars). In the fol- 
lowing sections, all parameter values that are functions of 
the source type S are given in Table 2 of W92. The fol- 
lowing sections give an overview of the distributions used; 
the reader is referred to W92 for more details. We note at 
the outset that the parameters characterizing these com- 
ponents, and in some cases even their existence (such as 
the molecular ring), are uncertain. In the discussion that 
follows, we will explore the effect of changing some of these 
components. 



3.1.1. The exponential disk 

The exponential disk has a density distribution given in 
cylindrical polar coordinates by 



p{R, z, S) = pd(S) exp 



R — i?i 



1*1 



h z (S) 



where h z {S) is the disk scale height, which depends on the 
spectral class S, Ro — 8.5 kpc is the distance from the Sun 
to the Galactic center, and h = 3.5 kpc is the radial scale- 
length. The exponential disk is truncated at i? max = 15 kpc. 



Arm 


a 


^max 




extent 


width 






(kpc) 


(rad) 


(rad) 


(kpc) 


1 


4.25 


3.480 


0.000 


6.00 


0.75 


r 


4.25 


3.480 


3.141 


6.00 


0.75 


2 


4.89 


4.900 


2.525 


6.00 


0.75 


2' 


4.89 


4.900 


5.666 


6.00 


0.75 


L 


4.57 


8.100 


5.847 


0.55 


0.30 


L' 


4.57 


7.591 


5.847 


0.55 


0.30 



3.1.2. The bulge 

The bulge has a density distribution given by 
p(R, z, S) = p B {S) x- 1 - 8 exp (~x 3 ). 

where 

x 



y/R? + k'f z 2 



Ri 

In the W92 model, the bulge axis ratio k\ is set to 1.6, and 
the bulge radius R\ is set to 2 kpc. 



3.1.3. The spiral arms 

The original W92 model included four main loga rithmic 
spiral arms, and one spur in the Solar neighborhood. ICohenl 
subsequently replaced the single local spur by two 



narrower spurs on either side of the Sun. The shape of the 
arms is parameterized as 



9{R) = a log 



R 



Rn 



where a is the winding constant, i? m i n is the inner radius, 
and # m in is the angle at the inner radius. The arms go out 
to .Rmax, but are truncated above a certain length. The 
par ameters for th e arms (compiled from Table 2 of W92 
and lGohenlll994fl are given in Table [TJ The arms follow 
the same radial and vertical exponential distributions as 
the disk, and have a constant width. Thus, their density 
distribution is effectively given by 



p(R,z,S) = p A (S)exp 



R — i?o 



h z {S) 



inside the arms, and p(R,z, S) = outside. As mentioned 
previously, the stellar population in this component is dom- 
inated by young, massive stars. 

3.1.4. The stellar (molecular) ring 

The SKY model contains a stellar ring with a density given 
by 



p(R, z, S) = p r (S) exp 



-(R- RrY 



2(7? 





~ 1*1 1 


exp 





where R r — 6.75 kpc is the ring radius, and oy = 0.96 kpc 
is the ring width. Since this component was motivated by 
an annulus of molecular gas in the inner galaxy, the stellar 
populations in this component are also weighted to young 
massive stars, as is the case for the spiral arms. 
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3.2. Dust distribution 

The second main ingredient for the Galactic model is a dis- 
tribution of dust, which is initially taken from the extinc- 
tion distribution of W92. For the dust properties described 
in i j3 . 31 the extinction distribution from W92 corresponds 
to 



p(R,z,S) = Pdexp 



h d z d 



(1) 



where hd — 3.5 kpc, Zd = 100 pc, and pd — 10~ 25 g/cm 3 . 

Given the dumpiness of the molecular gas of the Galaxy, 
the smooth exponential distribution of dust is clearly an 
oversimplification, but serves as a useful starting point for 
the model presented here. 

3.3. Dust properties 



The dust model adopted is from iDraine "fell (12001 . 
The model consists of a mixture of carbonaceous and 
amorphous silicate gr a ins w ith a size distribution from 
IWeingartner fc Draind (|2001h that reproduces interstellar 
extinction curves with Ry=3. 1. The model uses the renor- 
malization relative to H from iDrain 1 (|2003h . The size dis- 
tribution includes a population of small grains in two log- 
normal components in the size distribution, to represent 
transiently heate d very small grains and PAH molecules. 
The model from IDraine fc Lil used here is the one with a 
ma ss fraction o f PAHs set to qpAH=4.6%. 

IDraine fc Lil computed emissivities for dust grains il- 
luminated by the interstellar radiation field in the solar 
neighborhood from iMathis et all (|l983h scaled by a factor 
- denote d U - which varied between 0.5 and 10 7 . As de- 
scribed in IDraine fc Lil (|2007l ). the ionization fraction of the 
PAHs is a function of grain size. The emissivities and opac- 
ities used were identical to those publicly availabk[3 but 
were separated into three size regimes (B. Draine, private 
communication): ultra small grains (a < 20 A), very small 
grains (VSGs; 20A< a < 200A), and big grains (a > 200 A). 
These three regimes contain 5.86%, 13.51%, and 80.63% of 
the mass respectively in the dust model. Since the ultra- 
small grains are dominated by PAH molecules, we will refer 
to them as PAHs in the remainder of this paper (but we 
note that they are not strictly equivalent). 

3.4. Radiative transfer model 
3.4.1. Model set-up 

The radiative transfer models presented in this paper were 
computed with the Monte-Car lo three-dimens ional radia- 
tive transfer code Hyperioi\0 (Ro bitaillell201lD . This code 
can read in any arbitrary dust geometry, with multiple 
dust components, and multiple sources, making it perfectly 
suited to this project. It also allows images to be produced 
relative to an observer inside the density grid, which we 
make use of in this paper. The code implicitly solves the 
full equation of radiative transfer, including optical depth 
effects, scattering, polarization, and emission under the as- 
sumption of radiative equilibrium. 



The stellar and dust densities were discretized onto a 
cylindrical polar grid with 200 radial, 50 vertical, and 100 
azimuthal cells. The radial cells were distributed linearly 
between and 15 kpc. The vertical cells were logarithmi- 
cally spaced from 1 pc to 3 kpc above and below the mid- 
plane, with the region inside ±lpc divided into two cells. 
The azimuthal cells were uniformly distributed from to 
2tt. 

Rather than treating sources individually, which would 
be unfeasible computationally, 'diffuse' sources of emission 
were used to represent populations of stars: for each spec- 
tral class, each cell in the grid was given a probability for 
emission, and any photon emitted from that source was 
given the spectrum corresponding to the spectral class. 
Each diffuse source was given a total luminosity derived 
from the total number of 'real' sources. 

By default, the calibrated spectra used for each spectral 
class were derived from the absolute magnitudes given in 
Table 2 of W92 for BVJHK and 12 ^m & 25 pm. However, 
for m a in-seq uence stars, giant, and supergiants, we used 
lAllenl (|1973l ) to transform the spectral types into effec- 
tive temperatures and surface gravities, which we subse- 
quently used to int e rpola te stellar photosphere models from 
ICastelli &; Kuruczl (|2004j ). These model spectra were then 
scaled to the absolute magnitudes given in Table 2 of W92. 
The motivation for doing this is that the B magnitudes 
given in the SKY model are not sufficient to properly char- 
acterize the UV portion of the spectrum, which is crucial 
for the heating of the PAHs and VSGs. For the 'Young 
OB' source type in W9 2, which is meant t o represent OB 
associations, we fit the Caste lh~fc Kuruczl spectra to the 
magnitudes supplied by W92, and found that an effective 
temperature of 15,000 K provided the best fit to all wave- 
lengths. 

Each of the dust types (PAHs, VSGs, and big grains) 
was given the spatial distribution from Equation ([1} , or the 
modified spatial distribution defined in M4.2.21 scaled by the 
mass fraction of the dust types (c.f. t j3.3[) . 

3.4.2. Caveats 

The Hyperion radiative transfer code uses pre-computed 
emissivities for the grains for a range of interstellar radi- 
ation field strengths, rather than computing the transient 
heating of very small grains and PAH molecules exactly. As 
described in £13.31 the dust emissivities were pre-computed 
for a fixed illuminating spectral shape - but variable inten- 
sity U - so for regions in the Galaxy with different radiation 
field shapes, the shape of the emissivity will likely be dif- 
ferent. The primary effect of a different spectral shape is 
to change the amount of radiation absorbed by and there- 
fore heating the PAHs. Since PAHs have an opacity that 
is highest in the UV, harder radiation fields with a fixed 
intensity U will lead to higher excitation levels. This can 
be accounted for by quantifying the radiation field not by 
U (which ignores the wavelength dependence of the emis- 
sion) but by the power of the radiation field absorbed by 
the grains (per unit mass): 



.4 



1 http : //www . astro . pr inceton . edu/ draine/dust/irem . html where Jj sr f is the mean intensity of the radiation field. Thus, 

2 http://www.hyperion-rt.org two radiation fields with the same intensity U but different 
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Table 2. Variance between the models and the observations 





Original fH4.lt 


Two spiral 


arms (*I4.2.1|) 


Inner dust hole (U4.2.2D 


Increased PAHs (§4.2.31) 


Wavelength 


Longitude 


Latitude 


Longitude 


Latitude 


Longitude 


Latitude 


Longitude 


Latitude 


3.6 /im 


0.0133 


0.0056 


0.0071 


0.0019 


0.0057 


0.0008 


0.0064 


0.0047 


4.5 fim 


0.0248 


0.0097 


0.0179 


0.0062 


0.0155 


0.0009 


0.0113 


0.0002 


5.8 /im 


0.0756 


0.0246 


0.0235 


0.0046 


0.0255 


0.0093 


0.0112 


0.0011 


8.0 fim 


0.1147 


0.0542 


0.0374 


0.0142 


0.0444 


0.0242 


0.0187 


0.0016 


24 /im 


0.0643 


0.0095 


0.0331 


0.0041 


0.0323 


0.0027 


0.0406 


0.0074 


60 /im 


0.1191 


0.0560 


0.1640 


0.1027 


0.1163 


0.0331 


0.1061 


0.0251 


100 /xm 


0.0580 


0.0260 


0.0739 


0.0527 


0.0436 


0.0071 


0.0401 


0.0049 


Total 


0.0671 


0.0265 


0.0510 


0.0267 


0.0405 


0.0112 


0.0335 


0.0064 



spectral shapes - which would excite the PAHs by different 
amounts - would have correspondingly different values of A. 
Therefore, selecting the emissivities based on A is a better 
approximation than simply using U and takes into account 
the effects of spectral shape to first order. 

With this effect taken into account, the difference be- 
tween the emissivities due only to the change in spectral 
shape for a constant A is quan titatively smal l for our pur- 
poses. As shown in Figure 7b of|D ram 3 d20ll . for two very 
different radiation fields with the sam^l A, the difference in 
the PAH emissivities is less than a factor of 1.5. In this case 
the emissivities are systematically offset, and the colors do 
not change noticeably. Since at least some of the ISM will 
in fact have an interstellar radiation field close to that used 
to compute the dust emissivities, this effect may be smaller 
in reality. 

As shown in see Figure 7a of |D ram 1 (1201 lh . the dif- 
ference in the emissivities between fully ionized and fully 
neutral PAHs is more important. In particular, different 
levels of ionization can lead to differences in emissivities of 
factors of two, as well as significant changes in colors. 

We conclude that the assumption of the specific dust 
model used, which includes an assumption for the ioniza- 
tion level of the PAHs, is the main caveat in the radiative 
transfer models, while the method of choosing the emissivi- 
ties based on A, while approximate, appears to be adequate 
for our purposes. 

3.4.3. Model output 

Model images were computed for the survey area and col- 
lapsed into longitude and latitude profiles in the same way 
as the data. The latitude resolution of the profiles is 3', 
while the longitude resolution is 1°. The images were com- 
puted for 160 wavelength bins logarithmically spaced from 
3 to 140 microns, and were subsequently convolved with the 
transmission cu rves for IRAC, MIPS, and IRAS (following 
Appendix A of iRobitaille et al"1l2007t ). The convolution is 
especially important because the PAH features in the emis- 
sivities cause the flux to vary rapidly with wavelength in 
the mid-infrared, and picking a single wavelength instead 
of taking into account the proper transmission curve could 
result in the model fluxes being wrong by a factor of two 
or more. 

For all four models presented in Section [4] we have com- 
puted the variance between the data and the model in logio 

3 As described in Section 6 of iDraind (|201lD . the 20,000 K 
blackbody is adjusted in intensity to give the same power per H 
absorbed by dust, which is equivalent to saying with our defini- 
tion that the two radiation fields have the same A 



space to quantify the goodness of fit, and we list these in 
Table [5] We do not compute x 2 values, since the model 
is clearly misspecified - that is, the model is only a rough 
approximation of reality because it does not take into ac- 
count the small-scale features in the emission (which we 
are not attempting to reproduce) - so that absolute \ 2 va l~ 
ues and associated likelihoods would not be meaningful. 
We compute the variances separately for longitude and lat- 
itude profiles, and provide both the variances for individual 
wavelengths, and the overall values. 

4. Results 

4.1. Original model 

In this section, we show the results from the model as de- 
scribed in |J3]with no modifications. 

Figure Q] shows the average spectrum of the surface 
brightness inside -65° < I < 65° and -1° < b < 1° for 
both the observations and the model. The model under- 
estimates the surface brightness at IRAC wavelengths and 
overestimates it in the far infrared at IRAS wavelengths. 

In addition to the averaged surface brightness spectrum 
across the survey area, we also examine the spatial distri- 
bution of the emission. Figure [5] shows the observed and 
model longitude and latitude profiles from 3.6 to 100 /im. 
The emission in the model is of the right order of magni- 
tude, and reproduces a few important features of the ob- 
servations, including the scale height of the emission (long- 
wards of IRAC 5.8 /im) and the flattening of the IRAC 5.8 
and 8.0 /im emission inside \£\ < 30°, but there are also a 
number of issues: 

— the model surface brightness is too low by a factor of 
two or more for I < -30° for IRAC 3.6 /im to MIPS 
24 /im, and for I > 30° for IRAC 5.8 fim to MIPS 24 /im. 
Since £ = 30° corresponds to the tangency of the ring 
component for our model, this suggests that the model 
for the disk and/or spiral structure outside this radius 
are not adequate. 

— the model surface brightness is too high for \£\ < 30° 
for IRAS 60 /im and 100 /im, which suggests that our 
model contains too much dust inside R = 4 kpc. It has 
previously been suggested by numerous authors that a 
hole exists in the distri bution of gas within a few kpc of 
the Galactic center fc.f. lOrtiz fc Lepinell 19931 ) - presum- 
ably due to the bar - which would explain the excess 
flux predicted by the model here. 

— the IRAC 3.6 and 4.5 /im latitude surface brightness 
profiles are not as peaked as the data. This suggests 
that either the scale height of the stellar populations is 
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Wavelength (/im) Wavelength (jm) 

Fig. 1. Spectrum of the average surface brightness inside —65° < I < 65° and —1° < b < 1° for the initial model (left) and 
the final modified model (right). The thick colored solid lines represent different components of the model: supergiants 
(purple), AGB stars (blue), giants (teal), main sequence stars (green), other stellar types (dark gray), dust emission 
from PAHs (red), VSGs (orange), and big dust grains (yellow), and scattered light (black). The major components are 
labelled in the left panel. The shaded gray area shows the total flux. The open black circles show the model fluxes after 
convolution with the transmission curves, and the filled black circles show the observations. 



not small enough, or that dust extinction is too large in 
the mid-plane, since dust extinction is included in the 
model. The shape of the latitude profiles is dominated 
by the inner Galaxy (\£\ < 30°) since that is where most 
of the emission originates, so changes in this range will 
have the largest impact on the latitude profiles. 

4.2. Improved model 

In this section, we examine how the model presented in 
changes as we modify our initial assumptions in ways that 
have been motivated by previous investigations of Galactic 
structure. The modifications explored are changes to the 
parameters for the spiral arms (i j4.2.1l and Figure mod- 
ifications of the dust distribution interior to R = 4 kpc 
( £14.2.21 and Figure 2]), and modifications to the abundance 
of PAHs (i J4.2.3l and Figure [5]). We also explored the effect 
of concentrating the dust in the arms, but found that the 
simulations most consistent with the observations were the 
ones in which the dust is substantially broader than the 
stellar arms. 

4.2.1. Modified spiral structure 

Our modifications to the spiral structure model are moti- 
vated by the low model surface brightnesses for \£\ > 30° at 
IRAC and MIPS wavelengths relative to the observations 
Previous research has shown that when using stellar 
tracers in the near- or mid-infrared, there is only evidence 
for two of the spiral arms in the Galaxy (e.g. iDrimmel 
2000; IDrimmel fc Spergell [200lh . specifically the Scutum- 
Centaurus arm and the Perseus arm. On the other hand, 
surveys in the far-infrared and radio show that models with 
four arms are a better match. This suggests that the main 
potential of the Galaxy is in fact two armed, but that the 
non-linear response of the gas and dust flowing through 
the stellar potential is more complex, and includes arms, or 
large spurs, between the two major spiral arms. 



Table 3. Updated Spiral Arm Parameters 



Arm 


a 




$min 


extent 




width 






(kpc) 


(rad) 


(rad) 


(kpc) 


(kpc) 


1 


4.18 


3.800 


0.234 


6.00 


0.55 




1' 


4.18 


3.800 


3.376 


6.00 


0.55 




2 


4.19 


4.500 


5.425 


6.00 


0.55 




2' 


4.19 


4.500 


2.283 


6.00 


0.55 




L 


4.57 


8.100 


5.847 


0.55 




0.30 


L' 


4.57 


7.591 


5.847 


0.55 




0.30 



We therefore modify the model so that the Scutum- 
Centaurus and Perseus arms (2 and 2' in Table [T]) become 
the dominant arms by increasing the normalization factor 
of the O and B spectral class populations in these arms 
by a factor of two. Rather than completely eliminating the 
Norma and Sagittarius arms (1 and 1' in Table [T]), which 
does not provide a good fit to the dust-dominated bands, 
we eliminate all stars except those belonging to the O and 
B spectral classes, which we leave unchanged, under the 
assumption that stars are indeed forming in these secondary 
arms but do not stay there in the long term. 

The shape of these four arms is also changed to Gaussian 
profiles in the radial direction, with a — 550 pc. Finally, 
we modify the spiral arm parameters slightly to provide a 
better match to the shape of the longitude profiles. The 
updated spiral arm parameters are listed in Table [3l The 
local spurs are left unchanged. The resulting model lon- 
gitude and latitude profiles, shown in Figure [3] provide a 
better fit to the observations at IRAC and MIPS wave- 
lengths for \£\ > 30°, although we note that the IRAC 5.8 
and 8.0 /im fluxes are still slightly too low in that longitude 
range. This is reflected in the decrease in the variance be- 
tween the model and the data listed in Table [H for IRAC 
and MIPS wavelengths. 
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Fig. 2. The longitude and latitude surface brightness profiles in the range —65° < I < 65° and — 1° < b < 1° for 
the observations and for the initial model f M.ip . The thick black lines show the observed surface brightness, while the 
colored lines and the gray shaded area show the model surface brightness, with the same colors as used in Figure [T] The 
'traditional' spiral arm tangencies are indicated in the top panel. 
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Fig. 3. The longitude and latitude surface brightness profiles as for Figure [U but for the model including two major and 
two minor spiral arms instead of four main spiral arms ( H4.2.ip . 
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Fig. 4. The longitude and latitude surface brightness profiles as for Figure [21 with a deficit in the dust distribution in 
the central few kpc of the Galaxy ( N4.2.2[) . 
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Fig. 5. The longitude and latitude surface brightness profiles as for Figure BJ with 50% more PAHs f H4.2.3l) . 
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4.2.2. Modified dust distribution 

In order to provide a better fit to the IRAS 60 and 100 /im 
longitude profiles inside \£\ < 30°, we also experimented 
with including a hole in the dust within a few kpc of the 
Galactic center by modifying the radial dependence of the 
dust distribution from a simple exponential to a function 
with the following form: 




MoJ 



72a 2 



R ^ -^smooth 
R > -Rsmooth 



where /io and <tq are parameters, and /o and i? S mooth are 
defined as the normalization constant and transition radii 
for which the two functions transition smoothly, meaning 
that the functions and their derivatives are equal. These 
values are: 

al/h + no 



Rs 



and 



fo = 



3th 



exp 



smooth 



/h] 



exp 



smooth 



Mo J 



We found that values of /io = 4.5 kpc and oq = lkpc pro- 
vided a model that fit the IRAS profiles well, while main- 
taining the fit to the IRAC and MIPS profiles. A com- 
parison of the new density profile with the original one 
is shown in Figure The surface brightness profiles are 
shown in Figure [H and the corresponding average spectrum 
is shown in Figure [Q The MIPS 24 /im, IRAS 60 /an, and 
IRAC 100 /im now appear to fit the observed profiles well, 
while the IRAC 5.8 /im and 8.0 /im remain systematically 
low. 

As shown in Figure 2] and Table [5J the addition of 
the hole in the dust distribution also improves the fit of 
the IRAC 3.6 and 4.5 /im latitude profiles, which are now 
peaked in the same way as the observations, thanks to the 
reduced mid-plane extinction towards the Galactic center. 

4.2.3. Modified PAH abundance 

The fact that the systematic offset at IRAC wavelengths 
is seen at all longitudes suggests that it cannot be solved 
by changes in the geometry alone, but rather by a change 
in the large-scale properties of the dust or radiation field. 
For example, one solution might be to increase the abun- 
dance of PAHs relative to larger grains. To test this, we 
computed models with an increased dust density for the 
PAHs (strictly speaking the dust grains with a < 20A) , and 
found that a model with 50% more PAHs (corresponding to 
9pah ~ 6.7%) provides a good fit to the observation^. The 
resulting integrated spectrum is shown in Figure [T] and the 
profiles are shown in Figure [5] the fit to the IRAC 5.8 /im 
and 8.0 /im profiles is significantly improved (as seen in the 
decrease in the variance for these bands in Table [2]) , and 
the model now provides a reasonable fit at all wavelengths. 
This result suggests that - if the model for the stellar popu- 
lations is correct - the fraction of PAHs in the ISM may be 
higher than previously found, although we caution that due 



4 In reality, one should not increase the density of all grains 
with a < 20A, but instead, the qpAH fraction, which causes a 
bump in the grain size distribution, should be increased (since 
the 20 A threshold is arbitrary). 



a 
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-a 

> 
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Fig. 6. The relative radial dependence of the dust den- 
sity for the original model presented in £14 . 1 1 (gray) and the 
model with an inner dust hole from H4.2.2I (black) . 



to the large number of parameters involved, and given the 
caveats of the model (*j3A2j), this result is only tentative. 

We note that th e original value of (/pah = 4.6% found by 
iDraine fc Lil (pOfJl was also determined by using emission 
from th e Galactic pla ne, including data from the GLIMPSE 
survey. IDraine fc Lil implicitly assume that the emission 
and colors of the GLIMPSE emission are due to optically 
thin emission, but our model indicates that in the mid- 
plane, extinction can affect the IRAC fluxes by up to 50%. 
We ran a model without increasing the PAH abundance, 
instead computing the profiles as if the Galaxy was opti- 
cally thin, and found that the profiles and total flux were a 
better match to the observations, and fit almost as well as 
when increasing the PAH abund ance. Therefo re, it could be 
that our result is consistent with IDraine Sz Ll and that the 
true PAH fraction is higher than qpah = 4.6%, but appears 
to be lower due to the extinction of the IRAC wavelengths 
relative to the far-infrared. 

We also note that there may be other solutions to the 
systematic offset at IRAC wavelengths. For example, the 
clumpy nature of the ISM may affect the average energy 
input of the various dust types in a systematic way, result- 
ing in different relative contributions of emission from the 
different dust types at different wavelengths, but we defer 
such an investigation to a future study. 



5. Analysis 

In this section, we use the improved model as described 
in i j4.21 but we note that the qualitative results in ^5.11 and 
tj5.2l are not significantly changed if we use the initial model 
instead. 



5.1. Flux and energy breakdown 

Figures [1] to [5] show the breakdown of the surface bright- 
ness into contributions from the different components of 
the models, split up into main-sequence stars, giants, AGB 
stars, supergiants, other spectral classes, the three dust 
size ranges, and scattered light. We summarize this infor- 
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Fig. 7. The fractional contribution of various stellar populations and dust grain sizes to the total flux in the survey area 
as a function of wavelength. The colors are as in Figure [T] Only the main contributing spectral types and dust types are 
labeled. 
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Fig. 8. The fractional contribution to dust heating as a function of individual spectral type. The colors are as in Figure[TJ 
and the darker vertical lines indicate separations between individual spectral types. Only the main contributing spectral 
types are labeled. 



mation into a more intuitive representation in Figure [7] 
Approximately half of the IRAC 3.6 and 4.5 /zm flux is from 
giants, with the second and third biggest contributors being 
PAHs and main-sequence stars respectively. On the other 
hand, the IRAC 5.8 and 8.0 /Ltm bands are strongly domi- 
nated by emission from PAHs. MIPS 24 fim is dominated by 
emission two thirds from PAHs (in the form of continuum 
emission rather than features) and one third from VSGs, 
while IRAS 60 and 100 /zm are dominated by thermal emis- 
sion from the larger dust grains. 

A similar analysis can be done to understand which stel- 
lar populations contribute the most to the energy injected 
into and reprocessed by the dust populations. Figure [5] 



shows the breakdown of this energy by individual spectral 
type (with the main contributing spectral types indicated). 
Heating of the dust grains is dominated by B-type stars and 
to a lesser extent, OB associations. Since O stars are rarer 
and have shorter lifetimes, the heating from the O stars is 
likely to be less uniform, since O stars are mostly found in 
clusters. In fact, the strong peaks in emission seen in the 
longitude profiles, which correspond to regions of massive 
star formation, are likely to be heated by a larger fraction 
of O-type stars. For the larger dust grains, which are less 
sensitive to UV radiative than the PAHs and VSGs, giants 
provide around a third of the heating. 
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Fig. 9. The ratio of stellar emission measured in the GLIMPSE point source Archives to the total stellar flux predicted 
by the final model. This suggests that unresolved stellar flux contributes a large component of the total stellar flux inside 
\£\ < 20°. 



5.2. Unresolved flux 

The model derived in 34.21 can be used to understand how 
much of the diffuse emission in the IRAC bands is due to 
unresolved stellar flux. In order to study this, we first com- 
puted the mean surface brightness of point sources in the 
GLIMPSE Archive source lists in 1° bins in longitude, and 
compared this to the model stellar flux (excluding the dust 
emission) . Figure [9] shows the ratio of stellar flux measured 
in the GLIMPSE Archives compared to the model flux. For 
\£\ > 20°, the agreement between the resolved stellar flux 
and the model is good, suggesting that there is little or no 
unresolved stellar flux at these longitudes. Inside \£\ < 20°, 
the observed fluxes are as low as half of the predicted flux 
in places, suggesting that in that longitude range, the un- 
resolved stellar flux could be as important a component to 
the total stellar flux as the resolved stellar flux. 



5.3. External viewpoint 

As well as computing the appearance of the model from the 
position of the Sun, we can also compute the appearance 
of the Galaxy from an external viewpoint. In Figure ITOl we 
show what the model would look like at IRAC wavelengths 
viewed along the North Galactic pole. The spiral arms are 



prominent in PAH emission, while the bulge is bright in 
IRAC 3.6 /mi and 4.5 /iin. The molecular ring also features 
prominently, although as noted previously, there is debate 
as to whether this is in fact a real component of our Galaxy, 
or whether it is simply forme d by an inward extensio n of 
the major spiral arms (see e.g. lDobbs fc Burkert|[2?5Tl . 

6. Summary 

A radiative transfer model of the Galaxy which uses the 
SKY model i n con junction with the dust properties from 
iDraine fc Lil (j2007) was developed to self-consistently cal- 
culate the heating of dust grains, and verify whether it 
is able to reproduce the observed surface brightness from 
3.6 /im to 100 /im. The main findings presented in this pa- 
per are the following: 

1. The initial model is able to roughly reproduce the order 
of magnitude of the observed surface brightnesses ob- 
served, but there are disagreements between the model 
and observations, notably for \£\ > 30° for the IRAC 
and MIPS bands, and inside \£\ < 30° for IRAS 60 /im 
and 100 /im. 

2. By modifying the model to incorporate two major stel- 
lar spiral arms and two secondary spiral arms with only 
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x [kpc] x [kpc] 

Fig. 10. A color composite view of the Galaxy model from an external viewpoint, viewed along the North Galactic Pole, 
for the initial model (left) and the final modified model (right). The Sun is located at (x, y) = (0, —8.5 kpc). The colors 
show IRAC 8.0/xm (red), IRAC 4.5 //m (green), and IRAC 3.6 fim (blue), all on a square-root intensity scale. 



young massive stars, as well as removing dust from the 
central few kpc of the Galaxy, we are able to signifi- 
cantly improve the quality of the fit. A slight systematic 
offset remains at IRAC 5.8 and 8.0 /zm. By increasing 
the abundance of PAHs by 50%, we are able to elimi- 
nate this systematic offset, though we caution that other 
effects may explain this offset. 

3. The flatness of the IRAC 5.8 and 8.0 /zm and MIPS 
24 fj,m emission is not directly due to a hole in the dust 
distribution, but is a consequence of the lack of strong 
UV sources within the inner few kpc of the Galactic 
center. The evidence for a deficit of dust is only appar- 
ent at longer wavelengths, in the IRAS 60 and 100 /xm 
bands. 

4. Since the stellar populations are not symmetrically dis- 
tributed, the heating of the dust is not uniform, and 
therefore the spiral arms appear in diffuse emission 
at infrared wavelengths without the need for a non- 
symmetrical dust distribution (i.e. the model does not 
require additional dust in the spiral arms). This is not 
to say that the distribution of dust is not enhanced in 
the spiral arms, but that there is no evidence from the 
data that it is. 

5. The overall flux at IRAC 3.6 and 4.5 /jm is dominated by 
giant stars, PAHs, and main-sequence stars, while the 
flux from IRAC 5.8 /im to IRAS 100 /im probes mostly 
dust, with the size of the dust grains probed increasing 
with wavelength. 

6. On large scales, transiently heated very small grains 
and PAH molecules are predominantly heated by B-type 
stars. The larger dust grains also have a significant com- 
ponent of heating from giant stars. 

7. The Spitzer /IRAC bands may contain as much unre- 
solved stellar flux as resolved stellar flux for \£\ < 20°. 

The models presented here, while simple, help us quan- 
tify the main factors that determine the large-scale Galactic 



mid-plane emission. We plan to carry out significant im- 
provement to the models in the future, including a more 
systematic exploration of parameter space, a more realistic 
clumpy distribution of dust, and an improved treatment of 
the inner Galaxy, specifically relating to the molecular ring 
and the Galactic bar. In addition, the modeling of obser- 
vations can be extended to spatial and wavelength regions 
outside that covered here. For instance, one could include 
observations for the outer Galaxy (GLIMPSE 360) or at 
longer wavelengths (Herschel HiGal). Once the data from 
the Wide-Field Infrared Survey Explorer (WISE) are fully 
released, it will even be possible to model the all-sky ob- 
servations from 3.5 /im to 100 /im when combined with the 
IRAS all-sky data. 
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